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ABSTRACT 

An  experimental  investigation  of  the  internal  ballistics 
of  solid  fuel  ramjets  was  conducted  in  order  to  determine 
the  regression  rate  of  the  fuel  as  a  function  of  chamber 
pressure,  inlet  air  temperature,  and  air  flux  rate,  and  to 
model  the  flow  in  solid  fuel  ramjets  which  use  sudden  expan- 
sion flameholders  at  the  inlet.   In  addition,  flame  stabili- 
zation  limits  were  investigated.   A  computer  solution  for 
the  non-reacting  flow  field  gave  results  in  good  agreement 
with  experiments. 

Solid  fuel  ramjets  have  an  average  regression  rate  of 
the  fuel  that  closely  follows  the  theoretical  expression 
derived  for  kinetically  controlled  hybrid  rocket  combustion. 
The  inlet  step-height-to-motor-diameter  ratio  is  the  domi- 
nant parameter  in  determining  flame  stability  limits. 
However,  inlet  velocity  is  also  an  important  parameter. 
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I.   INTRODUCTION 

The  solid  fuel  ramjet  can  be  classified  as  a  special 
case  of  the  hybrid  rocket.   That  is,  the  fuel  is  solid  and 
the  oxidizer  is  gaseous  (air).   The  solid  fuel  ramjet  pro- 
vides many  of  the  advantages  of  solid  fuel  rockets  without 
the  weight  penalty  of  carrying  its  own  oxidizer.   Ease  of 
handling,  storability,  weight  savings,  and  possibly  greatly 
increased  tactical  ranges  are  among  its  advantages. 

A  research  project  is  currently  being  conducted  at  the 
Naval  Postgraduate  School  to  model  the  internal  ballistics 
of  solid  fuel  ramjets.   Such  a  model  would  provide  the 
design  engineer  with  a  working  tool  that  would  permit 
performance  parameters  and  operating  characteristics  to  be 
predicted.   The  overall  goal  of  the  research  project  is  to 
develop  a  model  that  can  predict  the  regression  rate  of  the 
fuel  as  a  function  of  operating  and  configuration  variables 
Very  little  experimental  data  currently  exist  that  can  be 
used  to  determine  the  adequacy  of  a  model,  thus  the  model 
must  be  developed  in  steps  verifying  each  step  with  experi- 
mental data.   Non-burning  tests  to  determine  flow  patterns, 
non-burning  tests  with  wall  blowing  to  determine  the  effect 
of  wall  mass  addition  on  the  flow  characteristics,  and  hot 
firing  tests  with  various  inlet  geometries  are  required. 
This  thesis  initiates  the  overall  research  project 
objectives. 
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Information  is  provided  from  several  areas  of  study 
that  are  pertinent  to  the  internal  ballistics  of  solid  fuel 
ramjets.   Specific  areas  are  flov;  and  heat  transfer  in-  ducts 
with  a  backward  facing  step,  flame  stabilization  in  high 
speed  flows,  hybrid  rocket  combustion,  and  thermal  decompo- 
sition of  polymers. 

Abbott  and  Kline  [Ref.  1]  investigated  the  turbulent 
flow  over  backward  facing  steps  using  water  flow  through  a 
two-dimensional  channel.   It  was  found  that  the  recircula- 
tion zone  consisted  of  a  complex  pattern  involving  three 
distinct  regions.   Immediately  downstream  of  the  step  face 
there  is  a  three-dimensional  zone  "characterized  by  one  or 
more  vortexes  rotating  about  an  axis  normal  to  the  vertical 
wall."   This  zone  is  followed  by  an  almost  two-dimensional 
zone  of  recirculating  flow.   The  area  at  the  flow  reattach- 
ment point  is  characterized  by  a  fluctuating  three- 
dimensional  flow.   It  was  also  found  that  the  location  of 
the  reattachment  point  remained  fixed  as  Reynolds  number 
and/or  inlet  turbulence  intensity  varied,  and  that  increas- 
ing the  step  height  moved  the  reattachment  point  further 
downstream  from  the  step  face. 

Krall  and  Sparrow  [Ref.  2]  using  water  flow  through  an 
orifice  in  an  electrically  heated  circular  tube  found  that 
local  heat  transfer  coefficients  in  the  separated,  reattached, 
and  redevelopment  regions  are  several  times  larger  than  those 
for  fully  developed  turbulent  pipe  flow.   They  found  that 


the  maximum  heat  transfer  coefficient  occurred  in  the 
vicinity  of  the  point  of  flow  reattachment. 

At  low  Reynolds  numbers  based  on  downstream  pipe  . 
diameter  (10,000  to  50,000)  the  heat  transfer  near  the 
reattachment  point  varied  rapidly  with  Reynolds  number. 
However,  at  higher  Reynolds  numbers  (70,000  to  130,000) 
only  a  weak  dependence  was  noted.   They  also  found  that  for 
weaker  separations  (small  step  heights,  h)  the  reattachment 
point  was  well  defined  but  for  strong  separations  (large 
step  heights)  the  reattachment  point  was  spread  out.   This 
transition  from  a  single  well-defined  reattachment  point  to 
the  "spread-out"  configuration  occurred  between  h/D=0.l6 
and  h/D=0.25,  where  D  is  the  downstream  pipe  diameter.   The 
reattachment  point  was  located  between  1.25  and  2.5  pipe 
diameters  downstream  of  the  initial  separation.   Increased 
step  height  moved  the  reattachment  point  slightly  downstream, 
Reynolds  number  and  Prandtl  number  had  no  significant  effect 
on  the  reattachment  point  location. 

The  results  of  Krall  and  Sparrow  are  also  in  general 
agreement  with  the  results  of  Abbott  and  Kline.   For  small 
step  heights  (h/D=.17)  the  two-dimensional  and  axisymmetric 
reattachment  points  are  very  nearly  equal  (X  /D  from  1.0  to 
1.25).   However,  for  larger  step  heights,  the  reattachment 
zone  spreads  out  considerably  more  for  the  two-dimensional 
case  than  for  the  axisymmetric  case. 

The  ideal  hybrid  rocket  is  one  in  which  a  gaseous 
oxidizer  flows  over  a  solid  fuel.   In  a  hybrid  rocket,  gas 
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phase  combustion  in  the  form  of  a  diffusion  flame  takes 
place  in  the  turbulent  boundary  layer  formed  by  the  oxidizer 
flow  over  the  solid  surface  of  the  fuel.   Heat  is  trans- 
ported by  convection  and  radiation  to  the  solid  surface 
which  causes  decomposition  of  the  fuel.   The  convective  heat 
transport  is  complicated  by  the  presence  of  high  rates  of 
wall  mass  addition.   Using  a  heat  transfer  limited  model, 
regression  rates  of  the  fuel  are  generally  proportional  to 
G  ,  where  G  is  the  total  gas  flux  and  n  is  a  constant 
(typically  0.8). 

Experimental  data  obtained  at  low  chamber  pressure 
and/or  high  oxidizer  mass  flux  indicate  that  the  regression 
rate  of  hybrid  grains  becomes  less  dependent  on  G,  but  is 
directly  related  to  the  partial  pressure  of  the  oxidizer  in 
the  gas  phase  [Refs.  3  and  4].   Reduced  reaction  rates  at 
lower  pressures  broaden  the  flame  zone.   Schlieren  studies 
have  also  shown  that  the  flame  zone  broadens  at  low  chamber 
pressures  and  may  broaden  enough  to  undergo  transition  to 
a  turbulent  premixed  flame  [Ref.  6].   When  the  oxidizer  is 
diluted,  the  diffusion  flame  moves  further  away  from  the 
surface,  increases  in  thickness,  and  the  flame  temperature 
drops  [Refs.  3  and  53-.   In  addition,  surface  kinetics  may 
be  rate  controlling  under  these  operating  conditions 
[Ref.  7].   When  these  conditions  exist,  the  heat  transfer 
limited  model  does  not  adequately  describe  the  combustion 
process. 
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Hybrid  rockets  operating  at  high  chamber  pressures  are 
adequately  described  by  a  heat  transfer  limited  model. 
Muzzy  [Ref.  4]  has  reviewed  the  modifications  that  have  been 
made  to  the  heat  transfer  limited  model  to  account  for  gas 
phase  and  surface  kinetic  effects.   The  modified  models 
[Refs.  6  and  7]  predict  a  dependence  of  the  regression  rate 
on  pressure  as  well  as  mass  flow  rate.   For  a  system  where 
radiative  heat  transfer  can  be  neglected,  Muzzy  has  shown 
that  the  regression  rate  equation  in  the  kinetically  domi- 
nated region  of  operation  can  be  expressed  as: 

r  =  C  P°-5  Ga  Xb 

where  C  is  a  constant,  P  is  pressure,  G  is  total  mass  flux, 
X  is  the  longitudinal  coordinate,  a  is  a  constant  that 
varies  from  0.3  to  0.4,  and  b  is  a  constant  that  varies  from 
-0.1  to  -0.2.   Another  study  by  Kustov  and  Rybanin  [Ref.  18] 
has  also  treated  the  pressure  sensitive  regime  of  operation. 
In  a  solid  fuel  ramjet  the  reacting  flow  has  character- 
istics that  are  similar  to  those  discussed  above.   However, 
two  significant  differences  exist.   Solid  fuel  ramjets  would 
normally  operate  at  lower  chamber  pressures  than  hybrid 
rockets,  resulting  in  significantly  greater  port  velocities. 
In  addition,  the  oxidizer  in  this  case  is  air  which  would 
correspond  to  a  severely  diluted  oxidizer  in  hybrid  combus- 
tion.  These  differences  require  some  means  of  flame  stabi- 
lization.  A  series  of  tests  performed  at  the  United 
Technology  Center  [Ref.  8]  have  shown  that  one  means  of 
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sustaining  combustion  in  a  solid  fuel  ramjet  is  to  use  a 
sudden  expansion  inlet  as  is  used  in  dump  combustors. 
Similar  results  were  reported  by  McCarthy  [Ref.  91. 

In  a  high  velocity  gas  stream,  where  the  velocities  of 
the  flow  are  at  least  an  order  of  magnitude  greater  than  the 
turbulent  flame  speeds  of  the  combustion  mixtures,  flame 
stabilization  is  required.   Bluff  body  flame  stabilization 
has  been  used  successfully  in  turbojet  and  liquid  fuel 
ramjet  combustor  design  in  the  past.   Typical  flame  stabi- 
lizers have  been  V-gutters  and  rods,  or  the  dome  of  a  can- 
type  combustor  [Refs.  10  and  11].   In  addition,  dump  or 
sudden  expansion  combustors  have  been  studied  for  integral 
rocket/liquid  ramjet  designs  [Ref.  191. 

The  recirculation  zone  behind  the  bluff  body  or  step 
causes  intense  mixing  of  fresh  unburned  fuel  and  air  with 
hot  combustion  products.   The  hot  combustion  products  ignite 
the  fuel-air  mixture.   This  stabilized  flame  can  then 
propagate  throughout  the  combustor.   Due  to  the  intense 
turbulent  mixing  in  the  recirculation  zone,  combustion  may 
in  some  cases  approach  that  of  a  well-stirred  reactor. 

A  well-stirred  reactor  model  considers  the  composition 
and  all  of  the  thermodynamic  properties  to  be  uniform 
throughout  the  volume.   An  idealized  reactor  with  instanta- 
neous mixing  would  have  homogeneous  fuel,  air,  and  product 
concentrations,  and  constant  pressure  and  temperature 
[Refs.  12  and  13]. 
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Solid  fuel  ramjet  combustors  have,  depending  upon  the 
combustor  length  and  inlet  configuration,  two  separate  and 
distinct  combustion  zones,  each  of  which  appears  to  be 
controlled  by  a  different  mechanism  [Ref.  9]  (See  Figure  1). 
In  the  recirculation  zone  products  and  reactants  are  rapidly 
mixed,  and  in  the  limiting  case  a  well-stirred  reactor  may 
be  approached.   Downstream  of  the  flow  reattachment  point 
the  boundary  layer  develops  and  the  combustion  may  be  similar 
to  that  in  a  hybrid  rocket  in  which  kinetic  effects  are 
dominant  (low  pressure,  high  mass  flux,  severe  oxidizer 
dilution) . 

Thus,  one  possible  model  for  solid  fuel  ramjet  combus- 
tion would  be  to  treat  the  recirculation  region  as  a  well- 
stirred  reactor  and  the  downstream  region  as  a  kinetically 
dominated  hybrid.   The  larger  the  grain  length-to-diameter 
ratio,  the  more  dominant  the  downstream  region  becomes  in 
determining  the  average  regression  rate,  and  the  more  one 
could  expect  the  ramjet  to  behave  as  a  kinetically  controlled 
hybrid  combustor.   However,  for  small  L/D  ratios  and/or 
large  recirculation  zones,  the  combustion  in  the  recircula- 
tion region  may  become  important  in  determination  of  the 
average  regression  rate.   In  addition,  the  combustion  char- 
acteristics in  the  recirculation  region  determine  the 
flammability  limits  for  the  entire  combustor.   Thus,  the 
model  used  in  this  region  requires  considerable  attention. 
For  very  strong  recirculation  zones,  a  well-stirred  reactor 
model  may  be  adequate.   However,  for  weaker  recirculation 
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zones,  which  are  desirable  from  a  pressure  loss  standpoint, 
kinetic  rates  may  be  of  the  same  order  or  greater  than  the 
reactant  mixing  rates.   In  this  case,  a  well-stirred  reactor 
model  cannot  be  expected  to  be  applicable  and  a  more  detailed 
consideration  of  the  recirculating  flow  is  necessary.   Thus, 
the  work  of  Gosman,  et  al.  [Ref.  14]  and  Spalding,  et  al. 
[Ref.  15]  becomes  increasingly  important  and  may  provide 
the  basis  for  modeling  the  recirculation  region  in  solid 
fuel  ramjet  combustors. 

It  is  the  purpose  of  this  study  to  determine  the 
dependence  of  the  regression  rate  of  a  solid  fuel  ramjet 
on  chamber  pressure,  inlet  air  temperature,  and  mass  flux 
of  air,  and  to  model  the  flow  in  the  recirculation  region 
caused  by  the  sudden  expansion  inlet. 
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II.   METHOD  OF  INVESTIGATION 

Non-burning  flow  visualization  tests  were  conducted 
using  a  plexiglass  model  with  three  inlet  step  configura- 
tions.  Flow  reattachment  points  of  the  incoming  stream 
were  measured  for  various  step  Reynolds  number  flows. 

A  Fortran  language  computer  program,  Pistep  II,  devel- 
oped by  Spalding,  et  al.  [Refs.  14  and  15]  provided  a  finite 
difference  solution  of  the  two  dimensional  Navier-Stokes 
equations  governing  heat  and  mass  transfer  in  recirculating 
flows.   It  was  used  to  calculate  the  flow  reattachment  point 
and  the  flow  field  for  non-burning  flows.   The  results  were 
compared  with  the  flow  visualization  data  and  the  work  of 
previous  investigators.   The  boundary  conditions  on  the 
differential  equations  were  altered  to  simulate  blowing  at 
the  wall  in  order  to  study  the  effect  of  mass  injection 
from  the  wall  on  the  flow  field. 

Experimental  firings  were  conducted  using  polymethyl- 
methacrylate grains  burning  with  non-vitiated  air  to  provide 
regression  rate  data  as  a  function  of  inlet  air  temperature, 
chamber  pressure,  and  oxidizer  flux  rate.   These  tests  were 
made  to  provide  initial  data  for  characterizing  ramjet 
combustion  and  flame  blow-off  limits,  and  to  determine  if 
the  reacting  flow  behavior  was  similar  to  that  obtained 
analytically  and  experimentally  for  non-burning  flow. 
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III.   DESCRIPTION  OF  APPARATUS 

A.  COLD  FLOW  VISUALIZATION  MODEL 

The  flov;  visualization  model  (Figure  2)  consisted  of 
three  sections.   The  first  section  was  a  2%   x  2^  x  5-inch 
plexiglass  block  with  a  three-quarter  inch  diameter  hole 
bored  through  the  center.   It  received  air  from  high  pressure 
tanks  and  supported  a  perforated  disk  flow  straightener . 
The  step  section  was  made  of  stainless  steel  and  was  ten 
inches  in  length  to  further  dampen  turbulence  and  provide 
fully  developed  velocity  profiles.   Three  steps  were  used 
with  inside  diameters  of  0.246,  0.519,  and  0 . 71^  inches 
respectfully.   The  third  section  was  a  15%-inch  plexiglass 
tube  with  an  inside  diameter  of  1.469  inches. 

For  the  lower  flow  rates  a  sonic  choke  with  a  diameter 
of  0.08l  inches  was  used  to  measure  flow  rates.   For  higher 
flow  rates,  an  orifice  flow  meter  with  inputs  to  a  Visicorder 
was  used  to  measure  flow  rates.   The  flow  reattachment 
points  were  measured  visually  using  a  scale. 

B.  HOT  FIRING  APPARATUS 

A  schematic  of  the  solid  fuel  ramjet  apparatus  showing 
the  overall  system  is  presented  in  Figure  3« 

1.   Ramjet  Motor 

The  solid  fuel  ramjet  motor  consisted  of  four  main 
sections:  the  head  end  assembly,  the  step  insert  section, 
the  grain,  and  the  aft  closure  section  (Figure  4).   All  of 
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the  metal  parts  of  the  motor  except  for  the  bolts  were  made 
of  stainless  steel. 

The  head  end  assembly  was  6k   Inches  long  and  contained 
the  inlets  for  the  air,  the  oxygen-methane  ignition  mixture, 
the  nitrogen  purge,  and  the  cooling  air.   Two  Champion  Z7 
spark  plugs  were  used  for  ignition. 

The  step  insert  section  was  2-3/4  inches  long  and 
held  the  step  in  place.   Three  steps  were  used  during  this 
investigation  with  inside  diameters  of  0.25,  0.50,  and  0.75 
inches.   This  section  also  accomodated  the  forward  end  of 
the  PMM  grain. 

Polymethylmethacrylate  was  chosen  as  the  fuel  in 
this  study  because  it  has  been  widely  used  for  basic  research 
in  hybrid  rocket  combustion  and  in  studies  of  polymer 
degradation  in  the  past.   It  is  composed  of  long  chains  of 
a  repeating  monomer: 


H      CH3 


H      C00CH3 


Idealized  chemical  reaction  with  complete  combustion 
would  involve  a  stoichiometric  reaction  with  air  as  follows 

C5H8°2  +  6°2  +  22-56N2  +   ^H20  +  5C02  +  22.56N2 
The  average  molecular  weight  of  the  products  in  this  case 
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is  31.3  lbm/lbmole.   The  gas  constant  is  R/M  =  49.4 
ft-lbf/lbm-°R,  and  the  stoichiometric  air  fuel  ratio  is 
8.27  lbmair/lbmfuel.   The  equilibrium  ratio  of  specific  heats 
at  a  temperature  of  3500°R  is  y   =  1.24.   An  adiabatic  flame 
temperature  calculation  for  the  stoichiometric  reaction  using 
a  value  of  H^  of  92.8  kcal/mole  for  the  heat  of  formation  of 
PMM  at  25  °C  [Ref.  5]  gives  a  flame  temperature  (with  no 
dissociation)  of  4265  °R. 

The  polymethylmethacrylate  grains  were  cut  from 
a  slab  into  3-3/4  x  3-3/4  x  12  inch  blocks.   The  ends  were 
rounded  to  an  outside  diameter  of  3-1/2  inches  to  fit  into 
the  step  section  and  the  aft  closure  section,  and  were 
sealed  with  "0"  ring  seals.   The  initial  inside  port  diameter 
of  the  grains  was  1.5  inches.   Eight  grains  were  made  that 
were  six  inches  in  length,  but  the  majority  of  the  tests 
were  performed  using  grains  that  were  12  inches  in  length. 

The  aft  closure  section  contained  a  straight 
section  with  a  pressure  tap  and  a  converging-diverging 
nozzle.   Nine  nozzles  with  various  throat  areas  were  used 
to  control  the  chamber  pressure.   A  Wiancko  0-200  psig  pres- 
sure transducer  was  mounted  on  the  forward  part  of  the  aft 
closure  section  to  measure  the  chamber  pressure. 

The  grains  were  held  in  place  by  four  3/8  inch 
diameter  steel  rods  between  the  step  section  and  the  aft 
closure  section.   Mounting  flanges  on  the  head  end  assembly 
and  the  aft  closure  section  were  bolted  to  a  large  test 
stand. 
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2.  Ignition  System 

The  ignition  system  metered  oxygen  and  methane  at 
800  psig  from  two  high  pressure  tanks  through  needle  valves 
into  the  head  end  assembly  section.   Inlet  ports  on  the  top 
and  the  bottom  injected  methane,  and  two  inlet  ports 
injected  oxygen  from  either  side.   Two  inches  downstream 
from  the  oxygen-methane  ports  were  two  Champion  Z7  spark 
plugs  located  180  degrees  apart.   A  transformer  supplied 
power  to  a  Model  T  spark  coil  to  energize  the  spark  plugs. 

3.  Nitrogen  Purge  and  Cooling  Air  System 

The  nitrogen  purge  and  cooling  air  system  was 
connected  upstream  of  the  head  assembly  section.   Two  high 
pressure  tanks  fed  nitrogen  into  the  motor  at  the  end  of  a 
hot  firing  for  approximately  one  second  to  extinguish  the 
combustion  flame.   Cooling  air  from  a  low  pressure  compressor 
was  then  blown  through  the  motor. 

k.      Air  Flow  Control 

A  standard  ASME  orifice  flow  meter  was  used  to 
measure  the  flow  rate  of  the  air  into  the  motor.   A  manually 
operated  gate  valve  between  the  orifice  and  the  motor  was 
used  to  provide  the  desired  flow  rate  through  the  motor. 
Two  pneumatically  operated  Jamesbury  ball  valves  operating 
together  either  vented  the  air  to  the  atmosphere,  or  allowed 
the  air  to  pass  through  the  motor. 

5.   Data  Acquisition  System 

A  Colvin  0-35  psi  differential  pressure  transducer 
was  used  to  measure  the  pressure  drop  across  the  orifice 
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flow  meter.   A  Wiancko  0-200  psig  pressure  transducer  was 
used  to  measure  the  pressure  upstream  of  the  orifice,  and 
an  iron-constantan  thermocouple  was  used  to  measure  the 
temperature  at  the  orifice.   A  second  iron-constantan  thermo- 
couple was  used  to  measure  the  temperature  of  the  air  enter- 
ing the  motor.   A  Wiancko  0-200  psig  pressure  transducer 
measured  chamber  pressure  at  the  aft  closure  section  of  the 

motor. 

\ 
The  orifice  thermocouple  output  was  recorded  on  a 

0-600°F  strip  chart  recorder.   All  of  the  other  transducer 

outputs  were  recorded  on  a  Honeywell  Model  2106  Visicorder. 

A  five  cycle  per  second  timing  signal  and  an  ignition  pulse 

were  also  recorded  on  the  Visicorder. 

6.   Air  Feed  System 

A  Pennsylvania  air  compressor  provided  air  at  pres- 
sures up  to  150  psia  at  a  rate  of  up  to  700  scfm.   This  air 
was  fed  into  an  air  reservoir  and  a  Polytherm  air  heater. 
The  Polytherm  air  heater  was  capable  of  yielding  1.75  pounds 
per  second  airflow  at  150  psia  and  1000°F.   The  heater 
burned  natural  gas,  and  the  air  passed  through  a  heat 
exchanger,  thus  providing  non-vitiated  hot  air. 

A  temperature  controller  directed  two  mixing  valves 
that  produced  the  required  air  temperature  at  150  psia.   A 
manual  gate  valve  controlled  the  air  flow  to  the  test  cell. 
A  valve  in  the  test  cell  controlled  the  flow  of  air  to  the 
pipe  leading  to  the  orifice  flow  meter. 

Figures  5,  6,  7,  and  8  show  the  test  cell  and  the 
ramjet  apparatus. 
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IV.   EXPERIMENTAL  PROCEDURE 

A.  COLD  FLOW  TESTS 

The  cold  flow  tests  were  performed  at  the  rocket  lab 
facilities  of  the  Naval  Postgraduate  School.   Visualization 
was  obtained  by  injecting  a  mixture  of  water,  alcohol,  and 
food  coloring  at  the  base  of  the  step  inlet  through  a  small 
tube  with  a  hypodermic  needle.   The  mixture  was  broken  up 
into  very  small  droplets  by  the  turbulence  in  the  recircu- 
lation zone.   In  the  recirculation  zone  the  droplets  on  the 
wall  moved  towards  the  step.   At  the  flow  reattachment  point 
the  droplets  on  the  wall  were  stationary.   Downstream  of  the 
reattachment  point  the  droplets  on  the  wall  moved  away  from 
the  step.   Using  this  technique,  the  flow  reattachement 
point  locations  for  various  step  Reynolds  number  flows  using 
three  step  sizes  were  measured  visually  using  a  scale. 

Initial  visualization  tests  were  attempted  using  a 
tufted  model  but  they  proved  to  be  of  limited  value  due  to 
the  highly  turbulent  nature  of  the  flow. 

B.  HOT  FIRING  TESTS 

All  hot  firing  tests  were  performed  in  the  jet  engine 
test  cell  at  the  Naval  Postgraduate  School.   Hot  firing 
data  were  obtained  from  a  series  of  firings  at  three  air 
flow  rates,  three  chamber  pressures,  and  three  inlet  temper- 
atures, using  three  step  heights  and  two  grain  lengths.   The 
majority  of  the  hot  firings  were  performed  with  a  12  inch 
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grain  length  and  an  inlet  diameter  of  0.5  inches.   Table  I 
summarizes  the  nominal  test  conditions. 

TABLE  I 
NOMINAL  TEST  CONDITIONS 


Chamber  Pressure  (psia) 
Air  Flow  Rate  (lbm/sec) 
Inlet  Temperature  (°F) 
Inlet  Diameter  (in) 
Step  h/D 
Grain  Length  (in) 


50 

75 

100 

.15 

.25 

.40 

70 

200 

350 

.25 

.50 

.75 

.  417 

.333 

.25 

12 

6 

For  the  tests  in  which  the  Polytherm  air  heater  was  used, 
the  temperature  controller  was  set  and  the  temperature  would 
stabilize  at  the  ramjet  inlet  in  one  to  two  hours.   This 
delay  was  due  to  the  fact  that  the  flow  rates  to  the  motor 
were  quite  low,  at  most  0.5  lbm/sec.   Once  the  temperature 
stabilized  at  the  desired  temperature  at  the  inlet  to  the 
ramjet,  the  manual  gate  valve  downstream  of  the  orifice  was 
set  to  obtain  the  desired  air  flow  rate. 

Chamber  pressures  were  preselected  and  controlled  by 
varying  the  nozzle  throat  area  for  a  given  air  flow  rate. 
Chamber  pressures  could  normally  be  attained  within  ten  psi 
of  the  desired  pressure.   Calculated  combustion  temperatures 
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varied  significantly  with  air  flow  rate  and  step  size.   Also, 
at  higher  chamber  pressures  the  preset  flow  rate  would 
decrease  upon  ignition  because  the  higher  back  pressure  would 
unchoke  the  air  flow  control  valve. 

Ignition  normally  lasted  for  four  seconds.   After  three 
seconds  of  ignition,  the  air  was  directed  through  the  ramjet 
motor.   Ignition  was  continued  into  the  first  second  of  the 
run  to  insure  that  combustion  would  be  sustained.   A  series 
of  thirteen  tests  at  high  flow  rates  in  which  the  ignition 
flame  was  blown  out  were  used  to  determine  the  rate  of 
consumption  of  the  PMM  grains  during  the  oxygen-methane 
ignition.   These  data  were  used  to  correct  the  initial 
weight  of  the  fuel  used  in  the  regression  rate  calculation. 

Combustion  lasted  for  thirty  seconds.   At  the  end  of 
each  run  the  motor  was  extinguished  by  simultaneously  venting 
the  air  to  the  atmosphere  with  the  Jamesbury  ball  valves  and 
actuating  the  nitrogen  purge  system.   Low  pressure  air  was 
then  blown  through  the  motor  for  cooling.   The  only  section 
of  the  motor  that  heated  up  significantly  was  the  aft 
closure  section  and  the  nozzle. 
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V.   PISTEP  II  COMPUTER  PROGRAM 

Pistep  II  is  a  Fortran  computer  program  developed  by 
Gosman,  Pun,  Runchal,  Spalding,  and  Wolfshtein  [Ref.  14], 
It  provides  a  finite  difference  solution  of  a  steady 
incompressible  turbulent  flow  through  a  cylindrical  pipe 
with  a  sudden  step  enlargement.   Equations  [Ref.  15] 
governing  the  distribution  of  stream  function,  ip,  vorticity, 
co/r,  non-dimensional  temperature,  T,  turbulent  kinetic  energy, 
k,  and  turbulent  energy  dissipation  rate,  £,  are  cast  in 
elliptic  form  in  cylindrical  polar  coordinates  as  follows: 

-  f-[b.p  * * —  ]  +  rSA   =    0 

8r    <j>   3r  <J> 

where  $   represents  any  one  of  the  variables  i^,  to/r,  k,  e, 
or  T.   The  functions  a,,  b  ,  C  ,  and  S   are  presented  in 
Table  II. 

A  two  equation  turbulence  model  is  employed  in  the 
program.   The  two  variables  are  turbulent  kinetic  energy,  k, 
and  the  turbulent  energy  dissipation  rate,  e.   The  effective 
viscosity,  u  ff>,  is  related  to  k  and  e  by 


"err  =  Vk2/e 
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A  "dissipation  length  scale",  £,  can  be  used  in  an 
alternative  formulation: 

1/2 
u  ff  =  C  kx/  p£ 
Keff    u     K 

Equating  the  two  equations  for  the  effective  viscosity  gives 

l  -   k3/2/s 

Constants  that  appear  in  the  two  turbulence  equations 

are  c-,  ,  c?)  c  ,  a^.,  and  a  .   Values  used  in  the  computer 

program  together  with  their  Fortran  symbols  are  given  in 
Table  III. 

TABLE  II 
VARIABLE  FUNCTIONS  (Table  I  of  Ref.  15) 


*         a*        b*        °*  s* 


2       2 
w/r      r       r      u 


eff 


2 
i[>       0     1/pr      1         -o)/r 


1    yeff/ak    X         -(G-pe) 


2 


1    ^eff^e    1         -(Cj^Ge/k  -  c2pe  /k) 


T      1    yeff/oh    1 


G  is  a  generation  term: 

dV^    0  3V^   0  V     0  av^        8V     p 

g  =  W2^>2  +  ^2  +  <-f>2i  +  i*r  +  ^  } 
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TABLE  III 
TURBULENCE  CONSTANTS  (Table  II  of  Ref.  15) 


Constant     c-,      c~      c       a.  o 


Fortran 
Symbol 


CI       C2       CMU     PR(NK)    PR(NE) 


Value        1.45    2.0     0.09    1.0      1.3 

Boundary  conditions  are  specified  at  every  point  on  the 
boundary.   Wall  functions  that  are  derived  from  Couette  flow 
solutions  are  used  in  evaluating  boundary  values  for  co/r,  k, 
e,  and  T.   Other  boundary  conditions  are  given  in  Figure  9. 

The  computer  program  calculates  heat  and  momentum 
transfer  in  steady  two-dimensional  axisymmetric  recirculating 
flows.   Five  simultaneous  coupled  non-linear  partial  differen- 
tial equations  are  solved.   The  finite  difference  equations 
are  implicit  equations  that  are  coupled  and  non-linear.   The 
Gauss-Seidel  point  iterative  method  is  used  to  solve  the 
finite  difference  equations.   Provisions  are  made  for  both 
overrelaxation  and  underrelaxation  of  the  variables. 

References  14  and  15  contain  a  more  complete  discussion 
of  the  mathematical  formulation  of  the  problem,  development 
of  the  finite  difference  equations,  and  solution  algorithms. 
An  excellent  discussion  of  several  turbulence  models  can  be 
found  in  Reference  16. 
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The  boundary  conditions  on  the  pipe  wall  can  be  changed 
so  that  the  effect  of  mass  addition  through  the  wall  for 
non-reacting  flows  can  be  investigated  analytically.   A 
constant  radial  velocity  at  the  wall,  which  yields  a 
linearly  varying  stream  function  along  the  wall,  is  used  to 
simulate  wall  mass  injection  (See  Figure  10).   In  this 
initial  study,  the  main  flow  and  the  wall  mass  addition 
were  identical  in  chemical  composition  and  temperature  . 

In  Pistep  II,  wall  functions  based  on  Couette  flow  are 
used  to  specify  the  boundary  conditions  on  the  pipe  wall 
for  all  'of  the  variables  except  stream  function.   For  very 
low  blowing  rates,  which  are  characteristic  of  solid  fuel 
ramjets,  Couette  flow  solutions  may  yield  realistic  results, 
However,  for  higher  mass  injection  rates,  Couette  flow 
solutions  may  not  adequately  represent  the  boundary 
conditions. 

The  centimeter-gram-second  system  of  units  is  used  in 
the  Pistep  II  computer  program. 
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VI.   RESULTS  AND  DISCUSSION 

A.   COLD  FLOW  EXPERIMENTS 

A  series  of  cold  flow  visualization  tests  were  conducted 

to  determine  the  reattachment  point  location  for  flow  over  a 

rearward  facing  step  at  the  inlet  of  a  cylindrical  tube. 

Defining  a  step  Reynolds  number  as: 

Re  ,    =  E where 

step      y 

p...  fluid  density 

V. . .  average  velocity  in  the  port  of  the  tube 

h.  .  .  step  height 

u.  . .  fluid  viscosity 

it  was  found  that  the  reattachment  point  location  remained 
constant  as  the  step  Reynolds  number  was  changed  for  a  given 
step  height.   The  reattachment  point  moved  further  downstream 
with  increasing  step  height.   The  tube  inside  diameter  was 
used  to  non-dimensionalize  the  reattachment  point  location 
from  the  face  of  the  step.   A  linear  least-square  curve  fit 
of  hj   data  points  was  used  to  obtain  the  following  equation 
(See  Figure  11). 


x 

-£  =  7.15  (h/D)  +  0.35 


Figure  12  compares  the  data  obtained  in  this  investigation 
with  the  data  of  Abbott  and  Kline  [Ref.  1]  and  Krall  and 
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Sparrow  [Ref.  2].   The  un symmetric  reattachment  points 
obtained  by  Abbott  and  Kline  are  shown.   The  flow  reattached 
at  different  locations  on  either  side  of  the  two-dimensional 
channel.   The  location  of  the  maximum  heat  transfer  coeffi- 
cient found  by  Krall  and  Sparrow  for  water  flow  through  an 
orifice  in  a  cylindrical  tube  was  used  to  locate  the 
reattachment  point  data  shown  in  Figure  12. 

The  data  from  this  study  fell  between  that  of  Krall  and 

\ 

Sparrow  and  the  downstream  reattachment  point  found  by 
Abbott  and  Kline.   Air  was  used  in  this  study  whereas  water 
was  used  in  the  other  studies;  however,  each  study  found 
that  reattachment  point  location  was  not  affected  by  changing 
Reynolds  number  for  a  constant  step  height.   Krall  and 
Sparrow  also  found  that  changing  Prandtl  number  of  the  flow 
did  not  change  the  reattachment  point  location. 

The  spread  in  the  cold  flow  data  in  Fig.  11  is  probably 
due  to  larger  relative  errors  in  measurement  closer  to  the 
step  face  and  not  due  to  Reynolds  number  effects. 

At  higher  step  Reynolds  numbers  small  vortices  were 
observed  at  the  step  face  similar  to  the  observations  of 
Abbott  and  Kline  for  two-dimensional  channel  flow.   The 
reattachment  point  was  axially  symmetric  for  most  Reynolds 
number  flows.   However,  at  high  Reynolds  numbers  the  reattach- 
ment point  varied  around  the  circumference  of  the  tube.   Four 
distinct  lobes  could  be  seen  when  this  occurred  and  may  be 
associated  with  the  small  vortices  at  the  step  face.   It 
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was  also  observed  that  the  rotation  and  mixing  of  the  flow 
in  the  recirculation  zone  became  more  intense  as  the  height 
of  the  step  increased. 

B.   PISTEP  II  RESULTS 

The  Pistep  II  computer  program  was  run  for  four  step 
configurations  and  five  inlet  velocities.   Table  IV  summarizes 
the  reattachment  point  results.   A  linear  least-square  curve 
fit  is  shown  in  Figure  11  and  Figure  12  for  comparison  with 
the  cold  flow  data  and  the  data  from  other  investigations. 
From  these  figures  it  is  seen  that  Pistep  II  gives  a  depen- 
dence of  X  /D  on  h/D  which  is  in  general  agreement  with 
experimental  results.   With  the  current  turbulence  model  and 
wall  functions,  the  results  agree  closely  with  those  of 
Krall  and  Sparrow  and  fall  somewhat  below  the  data  obtained 
in  this  investigation. 

A  21  x  21  grid  was  used  for  the  h/D  of  .25  and  a  grid 
of  20  x  20  was  used  for  the  other  values  of  h/D.   The 
computing  time  required  for  each  run  was  on  the  order  of  20 
to  30  minutes.   Underrelaxation  of  the  variables  was  required 
to  ensure  numerical  stability.   Even  with  underrelaxation, 
the  case  with  h/D  of  0.^21  was  numerically  unstable  for 
most  of  the  inlet  velocities  tried.   However,  this  case 
converged  the  fastest  when  numerically  stable. 

At  each  of  the  grid  points,  Pistep  II  calculated  the 
distribution  of  the  following  parameters:   vorticity,  stream 
function,  turbulent  kinetic  energy,  turbulence  dissipation 
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TABLE  IV 
RESULTS  OF  PISTEP  II  CALCULATIONS 


h/D 

X  /D 
r 

Re     * 
step 

Niter** 

Vin(cm/sec) 

.250 

1.735 

2446 

900 

1524 

.250 

1.696 

12232 

700 

7620 

.250 

1.761 

24465 

900 

15240 

.250 

1.702 

36698 

700 

22860 

.289 

1.979 

2009 

800 

1524 

.289 

2.008 

20089 

800 

15240 

.289 

2.004 

31033 

800 

22860 

.342 

2.152 

1335 

700 

1524 

.3^2 

2.199 

6677 

800 

7620 

.342 

2.161 

13354 

700 

15240 

•  342 

2.215 

20032 

800 

22860 

.421 

2.705 

4109 

523 

15240 

.421 

2.708 

6163 

521 

22860 

*  Based  upon  a  fixed  laminar  viscosity  of  I.78  x  10 
gm/cm. sec 

**  Niter  is  the  number  of  iterations 
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rate,  non-dimensional  temperature,  axial  velocity,  radial 
velocity,  effective  viscosity,  density,  and  turbulence 
length  scale.   Variations  of  shear  stress  and  heat  transfer 
coefficient  along  the  walls  were  also  calculated.   The 
location  of  the  flow  reattachment  point  was  determined  by 
solving  for  the  point  at  which  the  wall  shear  stress  was 
zero. 

C.   PISTEP  II  WITH  WALL  MASS  ADDITION 

The  boundary  conditions  were  altered  along  the  "north" 
wall  as  shown  in  Figure  10.   A  linearly  varying  stream 
function  (i.e.,  a  constant  radial  velocity  at  the  wall) 
was  used  to  simulate  injection  of  mass  uniformly  through  the 
wall.   For  this  initial  study,  the  chemical  composition  of 
the  mass  injected  through  the  wall  was  chosen  to  be  identical 
to  that  of  the  main  stream  coming  in  the  inlet.   FAR  is  a 
term  that  is  analogous  to  a  fuel  air  ratio,  and  represents 
the  fraction  of  the  total  incoming  inlet  stream  that  is 
injected  through  the  wall.   For  example,  FAR  =  0  gives  no 
wall  mass  injection.   FAR  =  0.5  gives  an  amount  of  wall  mass 
injection  equal  to  half  of  the  inlet  stream  so  that  the  amount 
leaving  at  the  exit  is  1.5  times  the  amount  entering  at 
the  inlet. 

FAR  was  varied  in  increments  of  0.1  for  the  case  in 
which  inlet  velocity  was  500  ft/sec  (VIN  =  15240  cm/sec) 
for  an  h/D  of  0.3^2.   Stream  function  output  was  punched 
onto  cards  which  were  used  as  input  for  a  separate  program 
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which  generated  contours  of  streamlines  using  the  Calcomp 
plotter  (see  Appendix  A).   Figures  13  through  21  show  the 
streamline  patterns  obtained  for  various  conditions.   It 
should  be  noted  that  the  streamline  patterns  are  shown  only 
for  the  upstream  half  of  the  ramjet  and  that  the  radial 
coordinate  has  been  expanded  more  than  the  axial  coordinate 
so  that  the  figures  are  not  to  scale.   This  permits  a  better 
visual  presentation  of  the  streamlines  in  the  recirculation 
zone.   Note  that  with  wall  mass  addition  (Figures  lH   through 
18),  a  streamline  that  originates  on  the  north  wall  in  the 
recirculation  region  escapes  from  the  recirculation  zone  and 
enters  the  main  stream  of  flow.   Figures  20,  135  and  19  show 
the  effect  of  increasing  step  size  (h/D)  while  maintaining 
inlet  velocity  constant  (mass  flow  rate  decreases).   Compari- 
son of  Figures  20  and  21  shows  the  effect  of  decreasing 
inlet  velocity  (or  inlet  mass  flow)  with  a  fixed  h/D  and 
comparison  of  Figures  13  and  21  shows  the  effect  of  decreasing 
step  size  with  a  fixed  inlet  mass  flow  rate. 

Table  V  summarizes  the  results  for  wall  mass  addition. 
The  reattachment  point  location,  based  on  the  point  where 
wall  shear  stress  is  zero,  moves  upstream  slightly  as  wall 
blowing  is  increased.   The  "wall  shear  stress"  is  calculated 
from  the  stream  function  difference  between  the  wall  node 
and  the  first  node  out  from  the  wall.   In  the  tube  between 
these  two  nodes  the  axial  velocity  is  at  least  an  order  of 
magnitude  greater  than  the  radial  velocity.   This  permitted 
the  calculation  of  the  "reattachment  point"  in  the  same 
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manner  as  for  no  blowing.   Also,  the  magnitude  of  the  radial 
velocity  at  the  wall  is  quite  small  compared  to  the  magnitude 
of  the  inlet  velocity.   It  should  be  emphasized  that  these 
results  were  obtained  using  the  same  boundary  conditions  on 
a)/r,  T,  e,  and  k  as  for  no  blowing.   Further  studies  are 
required  to  determine  whether  or  not  this  assumption  is 
realistic. 

TABLE  V 
PISTEP  II  WITH  WALL  MASS  ADDITION  (h/D  =  0.342) 


V.  (cm/sec)      FAR      x  /D      V  ^(cm/sec) 
in  r         wall 


15240 
15240 
15240 
15240 
15240 
15240 


0.0 

2.211 

0.0 

0.1 

2.191 

4.7 

0.2 

2.168 

9.5 

0.3 

2.145 

14.2 

0.4 

2.123 

19.0 

0.5 

2.101 

23.7 

Note:   900  iterations  for  each  of  the  above  runs 

D.   PMM/AIR  RAMJET 

Seventy-four  hot  firing  tests  were  conducted.   A  Fortran 
language  computer  program  was  written  to  reduce  the  data. 
A  computer  listing  together  with  the  output  is  included  in 
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Appendix  B.   The  first  three  pages  of  the  output  tabulate 
the  measured  data  for  each  test  in  which  combustion  was 
sustained.   This  is  followed  by  the  calculated  results. 
This  program  was  modified  to  reduce  the  data  for  the  tests 
in  which  combustion  was  not  sustained.   A  computer  listing 
and  output  for  the  non-sustaining  tests  is  included  in 
Appendix  C. 

The  expected  uncertainties  in  the  experimental  results 
was  calculated  using  the  error  analysis  method  of  Kline  and 
McClintock  [Ref.  17].   Using  the  notation  defined  in  the 
ramjet  data  reduction  computer  program,  the  uncertainties 
are  given  in  Table  VI. 

The  step  inlet  caused  non-uniform  regression  of  the  fuel 
grains.   There  was  less  regression  of  the  fuel  at  the  head 
end  in  the  recirculation  region,  and  nearly  uniform  regres- 
sion downstream  of  the  reattachment  point,  where  regression 
rate  was  observed  to  decrease  slightly  with  increasing 
distance  from  the  head  end.   The  inlet  effect  was  more 
pronounced  for  larger  values  of  h/D.   For  the  majority  of 
the  tests  there  was  a  significant  amount  of  carbon  deposited 
on  the  surface  of  the  fuel  at  the  head  end.   The  reattachment 
point  locations  were  measured  to  the  edge  of  the  carbon 
deposits  or  to  the  location  of  maximum  regression.   These 
two  locations  were  nearly  coincident  in  most  cases.   The 
regression  patterns  for  three  tests  are  shown  in  Figure  22. 

Regression  rates  were  calculated  using  two  methods. 
Based  upon  the  exit  diamter,  average  regression  rate  can 
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TABLE  VI 
ERROR  ANALYSIS 


Symbol 


%   Error 


Variable 


WDOT 

6.3 

TC 

13.8* 

RDE 

1.JJ 

RWT 

1.1 

WP 

.6 

PG 

8.6* 

PISP 

8.6* 

RESTEP 

5.8 

GAIR 

5.7 

UPORT 

15.0* 

SPC 

8.6* 

EQN 

11.0 

air  weight  flow  rate 

calculated  chamber 
temperature 

• 

r  based  on  exit  diameter 
r  based  on  weight  loss 
fuel  weight  flow  rate 
gross  thrust 
fuel  specific  impulse 
step  Reynolds  number 
air  flux  rate 
average  port  velocity 
specific  fuel  consumption 
empirical  r  equation 


*   The  uncertainty  in  the  ratio  of  specific  heats,  y>  and 
the  gas  constant,  R,  were  not  included  in  this  analysis. 
Complete  expansion  through  an  ideal  nozzle  was  assumed 
in  calculating  FG,  FISP,  and  SFC. 
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be  defined  by: 


•    df  "  di 

r  =   2  At —  (in/sec) 


where  d.  and  d«  are  the  initial  and  final  exit  diameters 
before  and  after  firing,  and  At  is  the  burn  time.   An 
average  regression  rate  based  on  weight  loss  of  the  fuel 
can  also  be  calculated: 


n 


4  (w,  -  w_)    P 

,  -  T  I  +  d .    -  d . 

r  =  2~At (in/sec) 


where  W.  and  Wf  are  the  initial  and  final  weights  of  the 
grain,  p  is  the  density  of  the  fuel,  and  L  is  the  grain 
length.   This  expression  for  the  regression  rate  gives  a 
better  average  value  when  there  is  non-uniform  regression 
along  the  length  of  the  grain,  and  was  used  throughout 
this  investigation. 

The  dependence  of  the  regression  rate  of  the  fuel  on 
chamber  pressure,  P  (psia) ,  inlet  air  temperature,  T  (°R), 
and  average  mass  flux  of  air,  G  (lbm/in  sec),  was  determined 
by  holding  two  of  these  parameters  at  fixed  values  and 
varying  the  third  for  several  tests.   It  was  found  that 
the  regression  rate  varied  as: 

•  .  cp0.51  T0.3*  „0.*1   (in/seo) 
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where  C  is  a  proportionality  constant.   For  a  step  h/D  of 

-4 
0.333,  C  =  2.3  x  10    and  for  a  step  h/D  of  0.416, 

-4 
C  =  3.4  x  10   .   A  plot  of  the  regression  rate  based  on 

weight  loss,  r   ,  versus  this  empirical  regression  rate 
equation  is  shown  in  Figure  23.   On  the  sixth  page  of  the 
ramjet  data  reduction  computer  program  output  in  Appendix  B, 
the  numerical  results  are  tabulated.   The  expected  uncer- 
tainty in  the  empirical  equation  was  calculated  to  be  11%. 

Muzzy  [Ref.  4]  has  shown  that  the  regression  rate 
equation  for  a  hybrid  rocket  operating  in  the  kinetically 
dominated  region  can  be  expressed  as: 


r  =  C  P0,5  Ga  Xb 


where  a  varies  from  0.3  to  0.4  and  b  varies  from  -0.1  to 
-0.2,  and  X  is  the  longitudinal  coordinate. 

Thus,  the  average  regression  rate  of  the  fuel  in  a  solid 

fuel  ramjet  is  very  similar  to  that  for  pressure  sensitive 

hybrid  rocket  combustion.   This  is  not  surprising  since  the 

reattachment  point  normally  occurred  at  less  than  one-fourth 
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of  the  grain  length.   Aft  of  the  reattachment  point  combus- 
tion should  approach  that  of  a  kinetically  controlled  hybrid 
combustor. 

The  maximum  regression  rate  of  the  fuel  occurred  near  the 
location  of  the  flow  reattachment  point.   This  is  due  to 
several  factors.   More  oxidizer  reaches  the  wall  from  the 
impingement  of  the  main  oxidizer  stream  at  the  reattachment 
point.   In  addition,  the  maximum  Nusselt  number  occurs  near 
the  reattachment  point  location. 

Data  from  36  of  the  tests  in  which  combustion  was  not 
sustained  are  tabulated  with  the  computer  program  in  Appen- 
dix C.   For  the  smallest  step  h/D  used  (0.25),  combustion 
was  sustained  only  while  the  ignition  system  was  operating. 
As  soon  as  ignition  was  terminated,  the  flame  blew  off. 
For  the  middle  step  (h/D  =  0.333)  the  flame  blew  off  as 
soon  as  the  air  flow  was  activated  for  many  of  the  tests 
at  higher  flow  rates.   Some  of  the  tests  in  which  combus- 
tion was  not  sustained  may  have  been  due  to  insufficient 
ignition  prior  to  activating  the  air  flow. 

Tests  performed  at  the  United  Technology  Center  [Ref.  8] 
have  indicated  that  inlet  configuration  is  more  important 
in  determining  when  combustion  will  be  sustained.   A  critical 
step  h/D  ratio  was  found  to  be  the  best  criterion.   For 
values  of  h/D  less  than  the  critical  value,  combustion  could 
not  be  sustained.   It  was  reported  that  for  one  particular 
all- hydrocarbon  fuel  the  critical  h/D  was  O.lU.   For 
polymethylmethacrylate  they  found  that  the  critical  h/D 
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was  0.31.   This  criterion  did  not  vary  with  step  Reynolds 
number. 

The  data  from  this  investigation  indicated  that  for 
high  inlet  velocities  with  h/D  greater  than  the  critical 
value  determined  by  UTC,  combustion  may  not  be  sustained. 
It  appears  that  not  only  is  h/D  important,  but  inlet 
velocity,  ignition  time,  and  inlet  temperature  may  also 
be  important.   At  identical  air  mass  flow  rates,  some  of 
the  tests  at  higher  inlet  temperatures  failed  to  sustain 
combustion.   This  is  probably  a  velocity  effect  since  at 
constant  mass  flow  rates,  increasing  temperature  increases 
the  velocity. 

The  Pistep  II  results,  although  for  non-burning  cases, 
should  indicate  trends  to  be  expected  in  the  burning  tests 
The  data  indicate  that  the  magnitude  of  the  velocity  at  a 
point  in  the  recirculation  zone  is  directly  proportional 
to  the  inlet  velocity.   Also  the  velocity  in  the  recircu- 
lation zone  V,  was  found  to  be  inversely  proportional  to 
h/D.   Thus  V  «  VIN  and  V  «  l/(h/D).   It  was  also  found 
that  -*r-  «  =-  but  did  not  vary  with  inlet  velocity. 

A  simplified  model  of  the  flame  stabilization  region 
is  one  in  which  the  recirculation  zone  is  considered  as  an 
ignition  zone  whose  length,  L,  corresponds  to  the  reattach- 
ment point  location.   At  the  flame  blowoff  condition,  the 
gases  traveling  through  the  recirculation  zone  have  an 
average  contact  time  equivalent  to  the  ignition  time  of 

the  gases:   t   =  t..   If  V  is  a  characteristic  velocity 
&        c    1 


in  the  recirculation  zone,  then  t  *   L/V.   Since 

3  c 


v   VIN 

V  « 


h/D 


and 


L  oc  x  /D  «  h/D 
r 


then  one  possible  expression  for  contact  time  is 


(h)2 

t    oc  IK- 

c    VIN 


The  step  h/D  has  a  more  dominant  effect  in  determining 
flammability  limits,  however  inlet  velocity  is  also  impor- 
tant.  This  treatment  is  oversimplified  and  was  discussed 
only  to  indicate  that  the  inlet  velocity  may  be  of  signifi- 
cance in  determining  blowoff  limits.   In  actuality,  ignition 
time  is  a  function  of  the  fuel-air  ratio  in  the  recirculation 
zone  and  this  is  in  turn  a  function  of  h/D.   Thus,  the  func- 
tional form  becomes  more  complex  as  one  considers  the 
details  of  the  recirculation  region. 

It  should  also  be  noted  that  the  simplified  model 
described  above  neglects  the  effects  of  turbulence  intensity 
on  flame  stabilization.   Examination  of  Pistep  II  output 
indicates  that  for  a  constant  inlet  mass  flow  rate,  increasing 
step  height  increases  both  the  turbulent  kinetic  energy  and 
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the  turbulence  dissipation  rate.   Also  for  a  fixed  inlet 
velocity,  increasing  step  height  decreases  the  recirculation 
zone  velocity,  the  turbulent  kinetic  energy  and  the  turbulence 
dissipation  rate. 

There  are  three  blow-off  situations  to  be  considered. 
The  first  one  is  a  steady-state  blowoff  limit  that  occurs 
when  the  mass  flow  rate  in  the  recirculation  zone  is  greater 
than  the  mass  consumption  rate  by  chemical  reaction.   The 
mass  flow  rate  through  the  zone  will  depend  upon  the  size 
of  the  zone  and  the  velocity  through  the  zone.   This  may  be 
important  in  throttling  applications  during  normal  combus- 
tion.  The  second  situation  is  one  in  which  there  are  igni- 
tion transient  blowoff  limits  during  starting.   Low  chamber 
pressure  during  ignition  gives  higher  inlet  velocities  which 
may  cause  the  flame  to  blowoff  even  though  combustion  would 
sustain  if  ignition  occurred  at  the  expected  operating 
pressure.   Thirdly,  it  is  possible  that  the  flow  at  the 
step  inlet  may  become  choked  at  high  mass  flow  rates  when 
the  chamber  pressure  is  low  during  ignition.   Subsequent 
expansion  could  modify  the  size  and  flow  characteristics 
of  the  recirculation  zone  and  cause  the  flame  to  blowoff. 
This  indicates  that  if  a  ramjet  operating  with  high  flow 
rates  could  be  started  at  a  higher  chamber  pressure, 
combustion  may  be  sustained. 

The  Pistep  II  computer  program  should  prove  to  be 
advantageous  in  studying  blowoff  limits.   When  chemical 
reaction  and  wall  mass  injection  are  included  in  the 
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program,  blowoff  velocities  in  the  recirculation  zone  could 
be  realistically  studied.   The  effect  of  h/D,  inlet  velocity 
(or  mass  flow  rate),  and  chemical  reaction  rate  could  be 
considered.   Blowoff  limits  have  been  calculated  in  this 
manner  for  premixed  gases  in  a  dump  combustor  [Ref.  15]. 

If  the  results  of  Abbott  and  Kline  [Ref.  1]  are  similar 
to  what  may  be  expected  in  an  axisymmetric  situation  with 
combustipn,  inlet  turbulence  would  not  be  expected  to  affect 
the  recirculation  zone.   However,  inlet  turbulence  may 
affect  the  regression  rate  of  the  fuel  downstream  of  the 
reattachment  point.   Further  tests  are  needed  to  clarify 
the  affects  of  inlet  turbulence  and  distortion  on  solid 
fuel  ramjet  combustion. 
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VII.   CONCLUSIONS 

The  cold  flow  visualization  data  were  in  general  agree- 
ment with  previous  studies.   Reynolds  number  had  no  measur- 
able effect  on  the  size  of  the  recirculation  zone  and  the 
reattachment  point  location.   The  reattachment  point  moved 
further  downstream  from  the  face  of  the  step  as  step  height 
increased. 

The  Pistep  II  computer  program  appears  to  predict  with 
acceptable  accuracy  the  reattachment  point  locations  for 
non-reacting  flows  with  no  wall  mass  addition.   Pistep  II, 
modified  for  wall  mass  addition,  indicated  that  the  size  of 
the  recirculation  zone  and  the  reattachment  point  location 
are  not  greatly  affected  by  wall  mass  addition.   The  location 
of  the  reattachment  point  moved  slightly  closer  to  the  step 
face  as  wall  blowing  was  increased.   Experiments  are  required 
to  verify  this  result,  and  the  assumptions  concerned  with 
the  boundary  conditions  used  in  the  computer  program  require 
further  investigation.   The  Pistep  II  computer  program  has 
proven  to  be  accurate  enough  to  warrant  further  development 
as  a  ballistics  program  for  the  recirculation  region. 

Solid  fuel  ramjets  which  use  PMM  as  the  fuel  have  an 

average  regression  rate  that  follows  the  expression: 

n  M    0  ?U    0  ill 
r  =  C  P   °   T   °   G  U*H±.   in  addition,  the  maximum 

regression  occurred  near  the  reattachment  point  location. 

Downstream  of  the  reattachment  point  the  regression  rate 
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decreased  slightly  with  increasing  distance  from  the  head 
end.   Thus,  solid  fuel  ramjets  have  an  average  regression 
rate  that  closely  follows  the  theoretical  expression  derived 
for  kinetically  controlled  hybrid  rocket  combustion. 

The  step  h/D  is  the  dominant  parameter  in  determining 
flame  stability  limits;  however,  inlet  velocity  is  also  an 
important  parameter. 

The  results  of  this  investigation  indicate  that  a 
reasonable  internal  ballistics  model  for  solid  fuel  ramjets 
may  be  obtained  by  using  the  Pistep  II  computer  program 
(modified  for  chemical  reactions  and  wall  mass  addition) 
together  with  the  regression  rate  equation  for  kinetically 
controlled  hybrid  rocket  combustion.   The  modified  Pistep 
II  program  should  be  used  to  analyze  the  recirculation 
region  (blowoff  limits,  etc.)  and  to  provide  the  input 
data  (chemical  composition,  etc.)  to  the  regression  rate 
equation  for  kinetically  controlled  hybrid  rocket  combus- 
tion.  The  latter . should  be  employed  downstream  of  the 
reattachment  point. 


V 
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VIII.   SUGGESTIONS  FOR  FUTURE  WORK 

This  work  was  an  initial  study  of  the  internal  ballis- 
tics of  solid  fuel  ramjets.   Additional  experimental  and 
analytical  studies  are  required  in  order  to  a)  obtain 
increased  combustion  efficiency,   b)  refine  the  internal 
ballistics  model  and   c)  study  the  effect  of  other  inlet  ge- 
ometries' and  flow  conditions  on  the  regression  rate  of  the 
fuel.   Specific  experimental  studies  needed  include  deter- 
mination of  a)  the  effect  of  wall  mass  addition  on  the 
recirculation  zone  for  a  step  inlet  using  inert  flows  (to 
compare  with  the  theoretical  predictions  using  Pistep  II), 
b)  the  effect  of  different  inlet  configurations  such  as 
swirl  vanes  and  plug  inlets  on  fuel  regression  rate  and 
flame  stability  limits,  and   c)  the  effect  of  different 
mixing  techniques  on  combustion  efficiency.   Analytical 
studies  required  are   a)  modifications  to  the  Pistep  II 
computer  program  to  include  chemical  reaction  with  heat  and 
mass  transfer  at  the  wall  and  b)  combination  of  the  recir- 
culation zone  model  (as  input)  with  the  kinetically  controlled 
hybrid  combustion  model  to  effect  a  usable  internal 
ballistics  model  for  solid  fuel  ramjets. 
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NORTH  WALL 


STEP  WALL 
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FIGURE  9.   BOUNDARY  CONDITIONS  FOR  PISTEP  II 
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